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Abstract 

The contribution of virtual s-channel Kaluza-Klein (KK) gravitons 
' to high energy scattering of the SM fields in the Randall-Sundrum 

■ (RS) model with two branes is studied. The small curvature option 

of the RS model is considered in which the KK gravitons are narrow 
low-mass spin-2 resonances. The analytical tree-level expression for 
Q^' a process-independent gravity part of the scattering amplitude is de- 

O . rived, accounting for nonzero graviton widths. It is shown that one 

cannot get a correct result, if a series of graviton resonances is replaced 
by a continuous mass distribution, in spite of the small graviton mass 
^ . splitting. Such a replacement appeared to be justified only in the 

^ I trans-Planckian energy region. 
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1 RS scenario with the small curvature 

The models with spacial extra dimensions (ED's) have pretensions to solving 
theoretical problems which have not yet received a satisfactory answer (such 

*Electronic address: alexandre.kisselev@mail.ihep.ru 



as hierarchy problem, proton hfetime, hierarchy of fermion masses and mixing 
angles, etc) and lead to a new phenomenology in the TeV energy region. The 
multidimensional gravity is strong, and the fundamental Planck scale can be 
related with the string scale. One manifestation of theories with ED's is the 
existence of Kaluza-Klein (KK) gravitons and their interactions with the SM 
fields. 

In the present paper, we consider one realization of the ED theory in a 
slice of the AdSs space-time with the following background warped metric: 

ds^ = e2«(-^-l?'l) dx" dx" + dy\ (1) 

where y = rO (— vr ^ ^ ^ vr), r being the "radius" of extra dimension, and 
ry^i, is the Minkowski metric. The parameter k defines a 5-dimensional scalar 
curvature of the AdSs space. 

We will be interested in the Randall-Sundrum (RS) model PP which has 
two 3D branes with equal and opposite tensions located at the point y = nr 
(called the TeV brane, or visible brane) and point |/ = (referred to as the 
Plank brane). If > 0, then the tension on the TeV brane is negative, 
whereas the tension on the Planck brane is positive. All the SM fields are 
constrained to the TeV brane, while the gravity propagates in five dimensions. 

The main goal of the paper is to estimate s-channel graviton contribution 
to the scattering of the brane fields in such a scheme. 

Let us note that the warp factor in the metric (P) is equal to 1 on the neg- 
ative tension (visible) brane, and a correct determination of particle masses 
on this brane is thus achieved [2]. By calculating the zero mode sector of the 
effective theory, one then obtains the "hierarchy relation" , 

Ml, = M (e2-- _ 1) , (2) 

K 

with Ms being a 5-dimensional Planck scale. 

From the point of view of an observer located on the TeV brane, there 
exists an infinite number of graviton KK excitations with masses 

mn = XnH, n = 1, 2 . . . , (3) 

where Xn are zeros of the Bessel function Ji{x), with 

x„ = 7r(n + i) +0(n-i). (4) 
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Note that all zeros of Ji{x)/x are simple ones, and that they are real positive 
numbers jH].^ 

The interaction Lagrangian on the TeV brane looks like (with the radion 
field omitted) 

1 1 °° 

Mp, ^ ^ ^ 

n=l 

Here T^'^ is the energy- momentum tensor of the matter on the brane, G^^J 
is the graviton field with the KK-number n, and 

A. = M, (6) 

is a physical scale on the TeV brane (here and in what follows, we neglect 
small corrections 0(e~'^'^'")). 

To get mi ~ 1 TeV, the parameters of the model are usually taken to be 
K ~ M5 ~ 1 TeV. Then one obtains a series of massive graviton resonances 
in the TeV region which interact rather strongly with the SM particles, since 

~ 1 TeV. 

In the present paper we will consider a different scenario which we call 
"small curvature option" jHE]: 

K < M5 ~ 1 TeV. (7) 

In such a scheme, the scale A^^ appears to be significantly larger than the 
gravity scale M5, as one can see from Eq. ®. Namely, we get A^ = 
100(M5/TeV)3/2(^100 MeV/ny/^ TeV. Nevertheless, for both real and vir- 
tual graviton production, a magnitude of a scattering amplitude is defined 
by the 5-dimensional Planck scale M5, not by A^^ or k, separately (see for- 
mulas in the next Section).^ 

Contrary to the case k ~ M5 ~ 1 TeV, there exists a series of very narrow 
low-mass spin-2 resonances with an almost continuous mass distribution. For 

^The minimal positive zero of the Bessel function Ju(z) is approximately equal to 

^Contrary to the KK gravitons, the radion will be hardly produced, since the radion 
couphng to SM fields is equal to l/(-\/3A^). The radion-Higgs mixing (if it exists, see 0) 
will be also suppressed. 
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such a case, the following inequalities were derived in Ref. [S]:^ 




(8) 



Notice, in order the hierarchy relation for the warped metric (j21) to be satis- 
fied, we have to put nr ^ 10. 

It is worth to underline that the AdSs space-time differs significantly from 
a 5-dimensional fiat space-time with one large ED even for very small k (i.e. 
for the small curvature). Indeed, let us consider the hierarchy relation for d 
compact ED's of the size Re'. 



where M4+rf is a gravity scale in the fiat space-time with d compact ED's. 
For d = 1, equation Q is a particular case of © in the limit 271 nr <^ 1. 
However, the condition 27rKr <^ 1 means that the ratio Ms/k should be 
unrealistically large: 



This inequality means, for example, that k -C 10"^^ eV, if M5 = 1 TeV. 

The present astrophysical bounds [Hj rule out the possibility d = 1 and 
significantly restrict the parameter space for d = 2, 3.^ The most stringent 
constraints come from neutron-star (NS) excess heat due to the trapped cloud 
of the KK gravitons surrounding the NS (for details, see the second paper 
in jH]). For instance, one gets > 4.4 ■ 10~^^ GeV (and, correspondingly, 
M4+1 > 1.6 • 10^ TeV) for d = 1. 

Fortunately, the above mentioned restrictions can not be directly applied 
to the AdS5 space-time, since they were derived in the soft radiation approx- 
imation, u; <^ T, where uj is the graviton energy, while T is the temperature 
of the nuclear medium (for instance, in the NS). Equation (jH)) means that 
a; ^ mi = 38(M5/TeV) MeV. On the other hand, a typical value of T is equal 
to 30 MeV. Thus, for M5 > 1 TeV the condition a; <C T is not satisfied, even 
for K = 10(M5/TeV) MeV. 

■^For the case k ~ AI5 ~ Mpi, analogous bounds look like 0.01 ^ k/Mp\ ^ 0.1 

one insists that M4+1 should be of order of few TeV, the case d = 1 is completely 
excluded, since Rc exceeds the size of the solar system. 



(9) 
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(10) 
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2 Virtual s-channel gravitons 

Let us study the scattering of two SM fields mediated by massive graviton 
exchanges in the s-channel, 



a a 



G^'^^^hh, (11) 



where a{h) = e~ , 7, q, g, etc. For instance, in hadron collisions virtual 
graviton effects could be seen in the processes pp 2 jets + X, pp 77 + X, 
and Drell-Yan process pp — ^ l^l^ + X. At linear colliders, the promising 
reactions are e~^e~ — >■ 77 and e+e^ — > //. In what follows, the invariant 
energy of the process, a/s, is assumed to be around 1 TeV.^ It means that 
we are working in the following region: 

> ~ M5 > K. (12) 

The matrix element of the process |TT| looks like 

M=AS. (13) 

The fist factor in Eq. p3|) contains the following contraction of tensors: 

A = T;, P^'^-P Tip = T;, T'^'' - ^ {TX iT% , (14) 

where P^''"^^ is a tensor part of the graviton propagator, while T^^"* is the 
energy-momentum tensor of the field a{h). 

2.1 Nonzero graviton widths 

We will concentrate on the second factor in Eq. (fT^ which is universal for all 
types of processes mediated by the s-channel exchanges of the KK gravitons. 
It is of the form: 

1 °° 1 

" n=l 

Here r„ denotes the total width of the graviton with the KK number n and 
mass rrtn- 



^The energy region ^ 1 TcV will be also briefly considered, see Eqs. H32(l - H35|l . 
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This sum was estimated in Ref. |3] in a zero width approximation (assum- 
ing r„ = for all n). It was shown that S{s) is purely imaginary in this limit 
(if no ultraviolet (UV) cutoff is imposed). Thus, there is no interference of 
ED effects with SM contributions to the same processes. 

The width of the massive graviton is indeed very small if its KK-number 
n is not too large jUj: 



rrin 



(16) 



where 77 ~ 0.09.^ However, the main contribution to the sum p5|l comes 
from the region n ~ \/s/k ^ 1. So, nonzero widths of the gravitons in the 
RS model with the small curvature should be taken into account. That is 
why, we will study a general case (r„ 7^ 0). For comparison, a particular 
case (all r„ = 0) will be also analyzed (see subsection l2l 
It is useful to present S{s) in the form: 



S{s) 
with 
Here 



1 



y I = - 

^ axi — hxl -\- c aia'^ — 



n=l 



Xt 



a = irjHi^, 6 = (kAtt)^, c = sA^ . 



xt 



(17) 



(18) 



a 



(19) 



and 



1 
2iri 



(20) 



are zeros of the quadratic equation az"^ — bz + c = 0. In the kinematical 
region (fT2|l . the parameter a |T9|l can be approximated as 



a 



K 



iri 

+ 1 



s 
M5 



(21) 



^To estimate r„, masses of the SM particles were neglected with respect to the graviton 
mass TO„ in]. This approximation is well justified since the sum in n is mainly given 



by the gravitons with the masses to„ ^ ^/s ~ f TeV. 
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with I (j| ^ 1. 

The sum in Eq. (|17|) can be calculated analytically by the use of the 
formula |2] 

where Zn,v {n = 1, 2 . . .) are zeros of the function z~'^Jiy{z). As a result, we 
obtain: 



ll-AlT] — 



1 _ 1 Up) 

a Ji(o-) p Ji{p) 



(23) 



Since rjs/Kl. <^ 1, we have the inequality: 

1 A, 



p| ^ ^ ^ > I a| , (24) 



and Eq. becomes 



with a given by Eq. ()21|1 (here and in what follows, small corrections like 
©(k/a/s) are omitted). The function S{s) has no singularities, as all zeros of 
Ji{z) are real, and Imcr 7^ at physical s. 

By using asymptotic behavior of the Bessel function jS] and formulas 
fjA.4|) in Appendix A, we obtain from (|25|) : 

1 sin 2A + z sinh 2e 

Sis) = „ ^ 15 — , 26 

^ ^ 4M|yi cosM + sinh^e ^ ^ 

where 

K 4' 2 VM5/ ^ ^ 

Formulas (j^^ and ()26|) are oitr main result. 

The following inequalities immediately result from (j^^ : 

- coth e ^ Im 5(s) ^ - tanh e (28) 
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Re5(s)| ^ 
ReSis] ^ 



1 + 2 sinh^e ' 



sinh2e ' 



(29) 
(30) 



where the notation S{s) = [2M^^/s] S{s) is introduced. Note that the upper 
bound for the ratio | ReiS(s)/ImiS(s)| decreases rapidly with energy, and it 
becomes as small as 0.08 at ^/s = 3M5. For comparison, this bound is equal 
to 0.85 at y/s = M5. The absolute value of Im4S(s) tends to 1 when s grows. 
For instance,_we find 0.98 ^ | Im5(s)| ^ 1.02 at y/s = 3M5. 

If ^/s ^ M5, the parameter e is numerically small, and we obtain that 



Sis) 



1 + y + O(e^) 



(31) 



at a/s = zqk, where zq is some zero of the Bessel function Ji{z). Thus, the 
value of S{s) at the point ^/s = zqk is actually defined by the graviton with 
the mass m„ = ^/s, while the relative contributions from other KK gravitons 
are suppressed at least by the factor rf / 12 ^ 7 ■ 10~^. 
At ultra-high energies, namely, at 

(32) 



(33) 



our sum can be approximated as follows: 



S(s) 



00 ^ 



n=l 



axl + c 



with the parameters a and c being defined above (|T8|) . The sum ()33|1 has the 
following asymptotics (see Appendix A): 



Sis) 



iB 



2\ 1/4 



where 



2V2 \t]s 

2\ 1/4 



TT 



(34) 



(35) 



Notice, the condition ^^frfs ^ A^^ means that we are working in the energy 
region ^ > 100 (M5/TeV)3/2(100 MeY / k)^'"^ TeV. 
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2.2 Continuous mass spectrum of the gravitons 

Since Xn — Tcn at large the mass splitting, AmKK — ttk, is very small with 
respect to the energy, and it seems reasonable to approximate a summation 
in n by integration over graviton mass rriKK- 

Let us calculate the real and imaginary parts of S{s) separately and 
compare the results of these calculations with Eq. (j26|) . We start from the 
following expressions: 



^ OO 



(36) 
(37) 



with a = s/h? and [3 = rfi^n/KT^Y. The usual way of calculating Eqs. (jHBj) . 
pTjl is to replace them by the integrals: 



Re Sis) 



1 1 



dz ■ 



dz 



1 - ^2)2 + 



(1 - ^2)2 + SZ^ 



(38) 



(39) 



where 



6 



/rjs 



A. 



< 1. 



(40) 



These integrals are estimated in Appendix B. The result of the calcula- 
tions is the following: 



Re5(s) 
ImiSfsl 



7]KS 



2k 



1 



^To be more correct, one should use ~ tt 

' an 



1 + 



8(irn) 



(41) 
(42) 

However, we can put 



nn at large values of n which are relevant for our calculations. 
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It follows from Eqs. (j4H) . (j42|) that the imaginary part dominates the real 
one. 

These expressions for the real and imaginary parts of S{s) do not agree 
with the "discrete mass" expression (j^Bj) . although they obey inequalities 
fl^ - ljnnj) . Remember that Eq. (pUj) was derived by the direct calculation of 
the input sum (jl5j) . However, the imaginary parts become practically the 
same in both cases in the trans-Planckian kinematical region, namely, at 
a/s > 3M5, as one can see from (j^HJ and ((121) • As for the real parts, they are 
small in comparison with the imaginary parts in this energy region.^ 

Can we approximate the discrete spectrum by the continuous mass distri- 
bution, if the mass splitting AmxK is "very small"? First of all, let us note 
that AmxK is dimensional and it should be compared with another dimen- 
sional quantity. Actually, we may regard a set of narrow graviton resonances 
to be a continuous mass spectrum (within some interval of n), if only 

AmKK < r„ (43) 

is satisfied. Let us stress, it is the inequality that allows one to replace a 
summation in KK number n by integration over graviton mass tukk-^ 

In our case, the relevant values of n, which give the leading contribution 
to the sought for quantity S{s), are n ~ y/s/ (ttk). Then we obtain from 
and (Uni): 

V^-^>nK, (44) 

or, equivalently, 

> 3M5 • (45) 
It is a common belief that in the fiat space-time with large ED's of the size 
Rc, the mass splitting is so small (AmxK = Rq^) ^^^^ the continuous mass 
approximation is undoubtedly valid. ^'^ Surprisingly, it is not a case. The 
reason is that the gravitons are extremely narrow resonances, r„ ~ m^/Mpy. 
Accounting for the hierarchy relation for d compact ED's Q, one finds from 
(I43p that only KK gravitons with unrealistically large masses, 

ml>Ml;--Mlll , (46) 



^See the comments and numerical estimates after Eqs. (I28ll - H30(l . 

^From the point of view of experimental measurements, the mass splitting must be 
compared with the experimental resolution Amj-cs- The spectrum looks continuous when 
AmKK < Amros, irrespective of Eq. (|43|l . 

i°For instance, R-^ « 130 eV, for d = 4 and M4+4 = 1 TeV. 
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are continuously distributed for d ^ 2. For all that, the widths of these 
gravitons remain relatively small, r„ > M4+(i (Mi+d/Mpi)^/'^. 
In particular, we get the conditions (for m„ ~ y/s): 



> M4+1 , 

v^> (MplM|+2)'/^ 



for d = 1, 
for d = 2, 
for ci > 1. 



(47) 



It is not surprising that the first of these inequalities is similar to Eq. (03)). 



2.3 Zero width approximation 

Now let us consider the limiting case of zero graviton widths (stable KK 
gravitons). The corresponding formulas can be obtained from the formulas 
derived above, if one formally takes the limit 77 — *■ 0, 77 > in them (remember 
that all r„ are proportional to rj). Let us introduce the notation 



S,(s) = Sis] 



r„=o 



Then we get from 



(48) 



P cot ( 

K 4 



71=0 



K 



TC[n + - 



(49) 



where V means the principal value. As one can see, neither So{s), nor 



Sis] 



I \/i=KZO 



dSH) depend on the large mass scale A^. 



This result ()49j) may be also obtained by direct calculations in the zero 
width approximation, if one uses the asymptotics of x„ at large values of 
n (jH) which are relevant in our case. Indeed, at y/s/n ^ 1 the quantity So{s) 
is approximated as 



Ms) 



E 

n=l 



1 



[■KK 



(50) 



and equation (jlHl) is then reproduced, as it is shown in Appendix C. 
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If we replace the summation in n by integration, we find that the imagi- 
nary part of iSo(s), 



5 n=l 



is of the form: 



oo 

Im5o(.) = jdm^iVS-m) = , (52) 



mo 



where mo = vr/t. This expression coincides with formula ()42|) derived in the 
previous subsection. 

The same procedure, when apphed to calculating Re5o(s), results in 

oo 

1 f 1 K 

Re>So(g) = -o V dm 5- = (53) 

mo 

The last term in (j53p is a particular case of expression (j4H) in the limit ^ 0. 
As one can see, iSo(s) is actually purely imaginary. Note that "continuous 
mass spectrum" formulas are in disagreement with the imaginary 

and real parts of the "discrete mass spectrum" expression taken at 77 — > 
(this discrepancy is discussed in the end of subsection 12. 2|) . 

Let us stress that a more rapid falloff of iS(s) at ^Jr\s ^ (see Eq. (jSH)) 
is completely lost in the zero width approximation. 



3 Conclusions and discussions 

In the present paper we have estimated the contribution of the virtual s- 
channel KK gravitons to the scattering of two SM fields. We have considered 
the small curvature option of the RS model with two branes (k <C Ms)- 
In such a scheme, the KK graviton spectrum is a series of rather narrow 
low-mass resonances. All the SM fields are confined to one of the branes. 

We have studied the case when both the colliding energy a/s and 5-di- 
mensional Planck scale M5 are equal to one or few TeV (and, consequently, 

^"'^The real part is suppressed with respect to the imaginary part by the factor k/ (2y/s). 
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^/s ~ M5 ^ k). Ultra-high energy region ^/s ^ M5 is also considered. By 
taking into account nonzero graviton widths (with n being the KK num- 
ber), we have derived tree-level analytical expression for S{s), the process- 
independent gravity part of the scattering amplitude (see our main formulas 
(123) and (Uni))- 

Then we have considered the case when a series of narrow low-mass gravi- 
ton resonances is replaced by a continuous mass spectrum. Both the real 
and imaginary parts of S{s) appeared to differ drastically from those derived 
without such a replacement. Thus, one has to conclude that an accurate 
estimation of the input sum in n is needed in order to get a correct result, 
epecially at < M5. The replacement of the summation in the KK num- 
ber by integration over graviton mass is well justified if the graviton mass 
splitting AmKK is less than r„. In its turn, this condition is actually satisfied 
only in the trans-Planckian region (when the invariant energy ^/s is several 
times larger than the 5-dimensional Planck mass). 

Zero width approximation (all r„ = 0) has been also studied. It is shown 
that some results can not be reproduced in this approximation, as one can 
see, for instance, from Eq. ()34j) which does not admit the limit r„ — >■ 0. 

Possible loop corrections to our results originating from the exchange 
of s-channel gravitons can be estimated in the zero width approximation 
as follows. Each new virtual graviton brings the coupling 1/A^, while a 
corresponding summation in n results in the additional dimensional factor 
1/k (after the substitution n —> mn/ini^)). Using dimensional arguments, 
we conclude, for example, that m-loop graviton contribution to a vector field 
scattering amphtude, A4^'^\s), should be proportional to 

2/i^3m-l 2/i^3m-l 

where M^ut is the UV cutoff in Feynman integrals, and m ^ 1. Thus, multi- 
loop effects may become dominating at Mcut ^ M^- In such a case, an 
effective operator analysis will be probably useful (see, for instance, Ref. j^). 
Similar arguments are also applied to the flat space-time with the large com- 
pact ED's, after replacements — > 1/M|i, 1/k ^ R'^. Note that the zero 
width approximation is very well justified for this case, since r„ ~ m^/Mpj 
(i.e. gravitons are extremely narrow spin-2 resonances, even for m„ ~ V^)- 

In a general case (nonzero r„), an additional nontrivial dependence on 
the scale appears in the RS scheme, since the graviton widths depend 



13 



on this mass scale (see Eq. Both simple dimensional arguments and 

formulas like Eq. (|54j) are no longer valid. 

Our formula ()25|) can be also applied to the scattering of the brane par- 
ticles, induced by t-channel graviton exchanges. Let 

Ml „ 

— > -t > /^^ (55) 
with t being 4-momentum transfer. Then we obtain from 

where lv{,^^ = exp(— zz/7r/2) J^iiz) is the modified Bessel function, and 

17] 



K 2 V M 



(57) 



Since l2{z)/Ii{z) ^ 1 at z > 1 (-7r/2 < arg^ < 37r/2), we find from ^ 
that 

in the kinematical region ()55|). Note that S{t) ()58|) is pure real and it co- 
incides with the imaginary part of S{s) derived in the zero width approxi- 
mation (j52|) up to the replacement s —t. In more general approach, one 
should sum KK-charged gravi-Reggeons, i.e. Regge trajectories a„(t) which 
are numerated by the KK number n ^21 H|- In such the amplitude 

has both real and imaginary parts. 

The relations of all cross sections necessary for studying effects induced 
by tree-level exchange of the KK gravitons with the quantities S{s) and S{t) 
can be found in the Appendix of Ref. 



Appendix A 

The sum (j33j) can be written as 

' ^ ' n=l 

^^See also Ref. [H] in which t-channel graviton contribution to the scattering amplitude 
was studied in a large curvature scenario of the RS model. 
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where a and c are defined by Eq. (fTH|). and 

C = e-'^l^ -i- =x-iy. (A.2) 

The sum in Eq. is of the form 

E°° 1 TT sinh -\/27rC + sin -\/27rC 1 /'a q\ 

n^ + C^ ~ 2V2C^ cosh V27rC- cos V27rC ~ 2^^ ' ^ ' ^ 

At large C, the main contribution to ()A.1|) comes from n ~ |C|. The 
disregard of term hx^^ in (fT7|) is justified if it is much less than c, that results 
in the inequality fevr^lCp <^ c, or, equivalently, ^Jrfs ^ A^^ ((221) • 

By using formulas 

sin(x — iy) = sin x cosh y — i cos x sinh 

cos(a; — iy) = cos x cosh y + i sin x sinh y, (A. 4) 

and 

sinh(a; — zy) = cos y sinh a: — z sin y cosh x, 

cosh(a; — iy) = cos ?/ cosh x + isiny sinh x, (A. 5) 

and taking into account that coshx ~ sinhx ^ coshy ~ sinh?/ is valid at 
X = cot(7r/8) y ~ 2.4?/ ^ 1, we obtain formula 



Appendix B 



At small ratio n/y/s, we get from Eq. that 

1 



ReiSnfsl 



^T^/sM| 



_ TTK 



where 



1= dz 



[l-z^)^ + 8z^ 



(B.l) 



(B.2) 



with 5 being defined by Eq. (j4(J|) . 
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It is convenient to divide our integral ()B.2|) into two parts: 

X = Ji + X2 . 

Here 



(B.3) 



du 



m(2-m) 



u{2 + u) 



_u2(2 - uf + 5(1 - uf u^{2 + n)2 + 6{1 + uY 



(B.4) 



and 



To = — du 



u{2 + u) 



u2(2 + m)2 + (5(1+^)8 ■ 
1 

Up to higher powers of 5, the integrals ()B.4|) and ()B.5|) are equal to 



(B.5) 



2 4 2 2v^ 

and we obtain formula (PT|) of the main text. 
As for the imaginary part, it is given by 



(B,6) 



where 



J = dz 



(B.7) 



(B.8) 



(1-^2)2 + ^^8 • 


As in the previous case, we divide the integral in ()B.8jl into two parts: 

J = Ji + J2, (B.9) 

with 

{l-uf (l + uY 



Ji 



du 



u^{2-uf + 5{l-uf u^{2 + uf + 5{l + uf 



(B.IO) 
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and 



1 



The leading parts of the integrals ()B.10|1 and (jB.lljl are equal to 



Thus, we get 



Zrj s 

that results in formula (jl^ in the main text. 

Appendix C 

Let us consider the sum 

M — n + f 

n=l ^ ' 

with 

M=^ + iO, v = ^ (C.2) 

(remember that ^/s ^ k). 

The infinite sum (ICLlj) can be found in Ref. [Til] : 

/C = ^ [vl/(t; - - + w)] - , (C.3) 

where "^{z) is the psi-function. Then one can apply the formula 

$(-^) = $(^) + i + 7rcot(7rz) , (C.4) 
z 

and use the asymptotic behavior of the function "^{z) at large z (| arg2;| < vr), 

M/(z)||^l^^^ = lnz-^ + 0(z-^) (C.5) 
17 



(both formulas are taken from PT). 

As a result, the asymptotics of /C looks like 



^ Ih»i = ^ - ^)] + ' (C.6) 



and we come to Eq. presented in the text. 
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